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■ ABSTRACT 

S3: 

We present intensive quasi-simultaneous X-ray and radio monitoring of the narrow 
^ ■ line Seyfert 1 galaxy NGC 4051, over a 16 month period in 2000-2001. The X-ray 

observations were made with the Proportional Counter Array on the Rossi Timing 
O ■ X-ray Explorer (RXTE) and radio observations were made at 8.4 and 4.8 GHz with 

jH ■ the Very Large Array (VLA). In the X-ray band NGC 4051 behaves very much like 

the analogue of a Galactic black hole binary (GBH) system in a 'soft-state'. In such 
systems, there has so far been no firm evidence for an active, radio-emitting jet like 
those found in 'hard state' GBHs. VLBI observations of NGC 4051 show three well 
, separated compact components almost in a line. This structure resembles the core 

^ ■ and outer hot spots seen in powerful, jet-dominated, extragalactic radio sources and, 

CT) | although no jet is visible in NGC 4051, suggests that a weak jet may exist. However 

CO . it has not previously been clear whether the nucleus is currently active in the radio 

■ band and whether there is any link between the radio and X-ray emission processes. 

Radio monitoring of the core of NGC 4051 is complicated by the presence of 

■ surrounding extended emission and by the changing array configurations of the VLA. 
' Only in the A configuration is the core reasonably resolved. We have carefully removed 
, the differing contaminations of the core by extended emission in the various arrays. 

The resulting lightcurve shows no sign of large amplitude variability (i.e. factor 50%) 
over the 16 month period and is consistent with being constant. Within the 6 A 
• t-h . configuration observations where we have greatest sensitivity we see marginal evidence 

■ for radio core variability of ~ 25% (~ 0.12mJy at 8.4GHz) on a 2-week timescale, 

correlated with X-ray variations. These percentage variations are similar to those of 
the Seyfert galaxy NGC 5548, which is 10 times brighter. Even if the radio variations 
in NGC 4051 are real, the percentage variability is much less than in the X-ray band. 
Within the B configuration observations, although sensitivity is somewhat reduced, 
there is no sign of correlated X-ray /radio variability. 

NGC 4051 is one decade lower in radio luminosity than the radio/X-ray /mass fun- 
damental plane for jet-dominated hard-state black holes, although the scatter around 
the plane is of the same order. The lack of radio variability commonly seen in hard 
state GBHs may be explained by orientation effects. Another possibility, consistent 
with the lack of radio variability, is that the radio emission arises from the X-ray 
corona although, in that case, the linear structure of the compact radio components 
is hard to explain. A combination of corona and jet may explain the observations. 
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1 INTRODUCTION 

Galactic black hole X-ray binary systems (BHBs) are often 
classified as being in one of a small number of 'states'. The 
most common states are the 'hard' and the 'soft' states. The 
hard state was originally defined as the state where the 2-10 
keV X-ray flux is low and the X-ray spectrum is hard, and 
the 'soft' state, where the 2-10 keV X-ray flux is high an d 
the X-ray spectrum is soft (Mc Clintock &; Remillardll2006l ). 
The hard state is dominated by a power- law component, and 
the 'soft' state is characterized by a strong thermal black- 
body component. The main physical parameter driving the 
states is the accretion rate wit h that rate being h i gher in 
the soft than in th e hard state ([Fender et al.l 120041 : iFenderl 
2010; Belloni 2010). Another major observational distinction 
between the hard and the soft states is that compact, power- 
ful radio-emitting jets are usually present in the hard states 
(jFendedl200H : IStirling etaD [2001) but have not yet been de- 
tected in the soft state. There have been detections of weak 
radio emission in s oft state systems (jBrocksopp et al J 12001 
ICorbel et all [2004) but this emission is usually thought to 
represent emission from relic plasma left over from a previ- 
ous hard-state jet and that no active jet is cu rrently being 
powered by the black hole ([Fender et al.ll2009h . 

NGC 4051 is a relatively nearby (15.21 Mpc, iRusselll 
2002) narrow line Seyfert 1 (NLS1) galaxy. It is one of the 
brightest and most variable Active Galactic Nuclei (AGN) in 
the X-ray sky and has bee n extensively observed by a num- 
ber of X-ray observatories (lM c Hardv et al|2004l ; IPonti et all 
120061 : iTerashima et al.ll2008l : iBreedt et aljboid ). One of the 
major results to emerge from these studies is that the X-ray 
variability characteristics of NGC 4051, as parameterised by 
the power spectral density (PSD), are very similar to those 
of the Galactic black hole X -ray binary (GBH) sy stem Cyg 
X-l when in the 'soft' state ( M c Hardy et al.ll2004h . A classi- 
fication of NGC 4051 as a soft state system is supported by 
the r elatively high accretion rate (~ 15% tHe) ([Woo &; Urrvl 
2002). If it were identical to a soft state GBH then an active 
radio-emitting jet would not be expected. Although the X- 
ray observations of other Seyfert galaxies are less extensive, 
the data are consiste nt with Seyfert galaxies being soft state 
systems (e.g. also see |M C Hardy et al.l ([20051 )). 

Radio variability has been detected previously in 
Seyfert galaxies but observations are not extensive. 
iNefffe de Bruvnl (|l983T ) detect variability of > 20% in the 
Seyf ert 2 galaxy Mkn 348 at 5GHz on ^yearly timescales 
and lWrobell (|2000h found variability at 8.4GHz in the Seyfert 
1 NGC 5548 of 33% ± 5% over 41 days and 52% ± 5% 
over 4 years. [Mundell et al.l (|2009h finds similar variations in 
most members of a sample of mainly Seyfert 2 galaxies and 
iFalcke et al.l ([2001 ) find the same in most members of sam- 
ples of radio- quiet and radio- intermediate quasars. However 
in none of the above cases we re there any quasi -simultaneous 
X-ray observations. Recently iBell et"aU {2010) have studied 
the X-ray and radio variability of the liner NGC 7213 and 
find a weakly significant correlation between those bands. 
Liners have much lower accretion rates than Seyfert galax- 
ies, being similar in accretion rate to the hard state GBHs. 
Liners are also much more radio loud than Seyfert galaxies, 
which is in agreement with the hard state scenario, and so 
an X-ray /radio correlation is expected. In soft-state GBHs, 
as the radio emission has been scarcely detected, we do 



not know what the exact relationship to the X-ray emis- 
sion might be. Correlated radio and X-ray observations of 
Seyfert galaxies may therefore provide us with important 
hints on the radio properties of the soft state systems 

The soft state Seyfert galaxy NGC 4051 has been known 
for some time to host a weak nuclear radio source of un- 
known origin (e.g. lUlvestad fc Wilson 1983 iKukula et al.l 
Il995l : IChristoDOulou et alJll997h . Based on the disappear- 



ance of the Hell 4686 line from therms spectrum at low 
X-ray flux levels, IPeterson et al.l (|2000h suggested that the 
inner part of accretion disc may change from being opti- 
cally thick to being advective at low accretion rates, lead- 
ing to the expectatio n of increased radio emission (e.g. 
iDi Matteo et al.1 Il999h . Therefore, over a period of 16 
months in 2000 and 2001, we carried out a joint radio (with 
the Very Large Array, VLA) and X-ray (with the Rossi X-ray 
Timing Explorer, RXTE) monitoring programme to search 
for an expected ant i- correlation between the radio and X-ray 
fluxes. 

Preliminary analysis of the resultant data did not, how- 
ever, reveal a strong radio/X-ray anticorrelation . Instead 
the datasets appeared to be positively correlated (|M c Hardvl 
2005). Correlated radio/X-ray variability is commonly seen 
in hard state, jet-dominated, GBHs, suggesting that the ra- 
dio emission in NGC 4051 may also come from a jet. We 
therefore undertook radio Very Long Baseline Interferom- 
etry (VLBI) observations both with the European VLBI 
Network (EVN) in 2003 and the Global VLBI network in 
2004. A triple structure consisting of a compact core with 
a barely resolved (~ few milliarcsec) component on each 
side of the core, separated by ^ 0.5 arcsec a nd almost in 
a straight line, was detected (|M°Hardvl l2005h . This struc- 
ture is similar to the core and hot spots of an FRII radio 
source and strongly suggests an underlying jet with the hot 
spots being the places where the jet hits some surrounding 
medi um. Similar EVN observ ations were later performed 
by iGiroletti &; Panessal (|2009h , revealing a broadly similar 
structure. A definitive analysis of the X-ray and radio vari- 
ability of the nucleus of NGC 4051 was, however, hindered by 
the fact that the VLA array configuration changed approx- 
imately every 4 months and thus the size of the synthesised 
beam changed, encompassing greater or lesser amounts of 
surrounding extended emission (see Fig[T]). In this paper we 
carefully model the radio emission of NGC 4051 on different 
spatial scales in order to remove the contribution from ex- 
tended emission to the flux from the nucleus, which we then 
compare to the X-ray flux. 

In Section 2 we present the radio observations and ra- 
dio analysis. In Section 3 we consider the radio variability 
as observed with the A and B array configurations, in which 
sensitivity to core radio flux variations is greatest, and com- 
pare it with quasi-simultaneous X-ray observations. We find, 
at best, only a marginal relationship. In Section 4 we discuss 
how to remove the contribution from extended radio emis- 
sion in the B, C and D array configurations. We find that, 
although the X-ray emission varies by factors of 10 or more 
over the 16 month observation period, there is no believable 
long term variation of the radio emission over the same pe- 
riod. In Section 5 we show simulations designed to determine 
the accuracy with which the extended contribution can be 
measured. In Section 6 we consider the implications of our 
result for radio emission models for radio quiet AGN. We 
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include discussion of jet models which underly the so-called 
'fundamental pla ne' relating black hole mass to radio an d X- 
ray luminosity (|Falcke et alJl2QQ4l: iMerloni et aHfe oOS) and 
we also discuss coronal models (Laor & Behar 2008). 



Table 1. Observational Details of VLA data at 4.8 and 8.4 GHz. 
For both frequencies the phase calibrator is J1219+484 and the 
flux calibrator 3C286. The phase calibrator was generally ob- 
served for £^5 min and the flux calibrator for « 4 min per run. 



2 RADIO OBSERVATIONS AND DATA 
REDUCTION 

NGC 4051 was observed by the VLA 29 times between 2000 
16th June and 2001 3rd September with a typical interval be- 
tween observations of 2 weeks. On each occasion the source 
was observed at both 8.4 and 4.8 GHz with a total time on 
source, in each band, of approximately 12 minutes. Six of 
our observations were made with the VLA in A configura- 
tion, 2 in BnA, 6 in B, 1 in BnC, 6 in C, 3 in DnC and 
5 in D. We observed in 2IF mode with 50MHz bandwidth 
at both frequencies. The phase calibrator was J 12 19+484 
and the flux calibrator was 3C286. For more information on 
these observations see Table [U 

The data were flagged and calibrated, using MPS, in 
the standard manner. Maps were made using imagr with 
clean fields set around two fainter neighbouring sources, one 
approximately 5.5 arcmin to the NE and another approxi- 
mately 3.5 arcmin to the south. After experimenting with 
different iterations of cleaning we standardised on 10,000 it- 
erations. Self calibration was not used on the final images. 
All early experimentation with self calibration revealed the 
source to be too faint, and resulting images were not im- 
proved by this method. 

Our aim here is to search for a possible relationship 
between the X-ray and radio fluxes. The X-rays vary very 
rapidly, and hence come from a very small region, tens or 
hundreds of light seconds across. If there is any relationship 
with the radio emission, we must therefore expect the radio 
emission to come from a region which, although possibly not 
as small as that of the X-ray emission region, is unlikely to be 
light years across. At the distance of NGC 4051 (~ 15Mpc), 
1 arcsec^ 70pc and hence we are concerned here with the 
core radio flux. However measurement of that flux accurately 
for a source with any extended structure is far from trivial 
and, even for an almost unresolved source, it is not easy, as 
we discuss further in Section |3J"1 

The standard AIPS task for measuring the peak and in- 
tegral flux densities is jmfit. This task fits a gaussian model 
to the values within a box drawn around the source and out- 
puts the integral flux density (Jy) values and peak flux den- 
sity (Jy/beam) values. Although we are interested in intrin- 
sic core flux, measurements of both peak and intergral flux 
densities have value for our study and so both are recorded 
and listed later. 

A box was placed around the nuclear core of NGC 4051 
in each map using the task COWINDOW. The box size was 
different for each array configuration and was chosen so it 
included all the nuclear flux. We kept the box size constant 
when calculating flux density values within the same array. 
jmfit outputs errors depending on the gaussian fitting to the 
nucleus of the source. These jmfit errors are then combined 
in quadrature with the standard 5% error on the flux density 
value to calculate the error on each of the 29 flux readings. 
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3 VARIABILITY OF THE NUCLEUS 

3.1 During period of A configuration observations 

3.1.1 8.4GHz Observations 

In Fig[T]we present images of NGC 4051 at 8.4GHz from rep- 
resentative individual observations made in each of the four 
VLA array configurations. We can see that, with the A con- 
figuration, the core is separated from the neighbouring east- 
ern and western components. Higher r esolution maps made 
from observations by MERLIN (e.g. IChristopoulou et al. 
1997h ) or the EVN (|Giroletti fe Panessal 120091 : lM c Hardv 
2005) show the same structure. Thus the flux of the central 
component, as seen in the A configuration observations, can 
be taken as reasonably representative of the nuclear flux. 
In all other array configurations measurements of the nu- 
clear flux are contaminated by extended emission, which re- 
stricts the conclusions which can be drawn from these data. 
We consider the problem of contamination in more detail in 
Sections 3.2 and 3.3 but here consider the specific concerns 
related to the A configuration observations. 

The integral flux density is the desired measurement 
for slightly extended, but otherwise isolated sources, such 
as one might find in deep surveys. However the peak flux is 
generally the best measure of the true core flux. If the source 
is compact and isolated but there are significant phase er- 
rors in the calibration which move some of the flux slightly 
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Figure 1. Representative maps made from single 8.4GHz observations in each of the four main VLA array configurations. The configu- 
ration (A, B, C or D) is listed in the top left corner of each map. The maps are all centred on the nucleus of NGC 4051. The rms noise 
level in each A, B, C and D map is 16, 2, 2.2 and 1.61/xJy/beam and the major axis of the restoring beam is 0.22", 0.78", 2.54" and 7.83" 
respectively . The contours are at rms x -2.8, 2.8, 4, 5.6, 8, 11, 16, 23, 32, 45, 64, 90, 127, 180, 254. 



to the side, then the integral flux density would be a bet- 
ter measure. (For a very brief discus sion of how one m ight 
measure core fluxes see, for example, iMullin et al.ll2008l and 
papers referred to therein.) However if the core is surrounded 
by extended emission then the integral flux density will be 
much more affected by small changes in beamshape which 
may include different amounts of extended emission. Our A 
configuration maps, although showing a core which is close 
to being unresolved (eg Fig [T]) consistently indicate an ex- 
tension, of approximately 0.1 arcsec, along the line of the 
3 main source components and so the problem of extended 
emission exists even for the A configuration observations. 

In the top panel of Fig [2] we show the 8.4GHz peak core 
flux densities with fixed beam (filled circles) and default 
beams (open hexagons) as derived using jmfit. For the non 
fixed beam maps the restoring beams were left free, to be 
determined only by the UV coverage and the use of robust 



in IMAGR. The size of the fitted gaussian component was also 
left free. As the observations were made within a relatively 
short period, at mostly the same hour angles, the restoring 
beams, in practice, were very similar, although not identical. 
The fitted gaussians were similar, although there was some 
scatter in width. 

In order to try and standardise further, we remade all 
the maps with a fixed restoring beam with parameters at 
the mean of the unconstrained restoring beam (i.e. for A 
configuration beamsize 0.24" x 0.19" and position angle of 
137 degrees) and robustness 0. 

The first point to note from Fig [2] is that any variation 
of the core 8.4GHz flux, if real, is very small, and certainly 
of much lower amplitude than that seen in the X-ray ob- 
servations (lower panels of Fig [2}. The quoted errors are 
the statistical error given by jmfit combined in quadrature 
with the 5% error usually assumed as the maximum likely 
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uncertainty on the flux density calibration. The calibration 
uncertainty therefore slightly dominates the total resultant 
error which is typically twice the statistical error, or slightly 
more. This total error is the correct error to use when com- 
paring the flux density from the same object in different 
observations. To determine the relative flux densities of dif- 
ferent objects within the same observation, however, one 
would not include the flux calibration uncertainty. 

We also note that minor calibration errors, or sidelobes 
from distant sources which have not been perfectly removed, 
can typically lead to flux measurement errors which are 
larger by factors of two, or maybe more, than the statisti- 
cal error. We can use observations of other structures in the 
maps to make a crude estimate of the real errors on the flux 
density of the central source. In our case components are vis- 
ible on either side of the core in the A configuration maps 
which we might assume to be constant if they arise from 
extended components related to NGC 4051. The (brighter) 
eastern component is definitely more extended than the core, 
and the (fainter) western component is probably the peak of 
a low surface brightness extended structure which extends 
to the south and so flux measurements of both will be a 
little more sensitive to the exact beamshape than will mea- 
surements of the core. Nonetheless, the default beamshape 
does not change greatly but variations in peak flux by twice 
the statistical error, and occasionally by more, are seen be- 
tween observations. We conclude that the error which we use 
throughout this paper, i.e. the statistical error from jmfit 
combined in quadrature with a 5 percent amplitude calibra- 
tion error , may be a slight underestimate of the true total 
error but probably not by more than 30 percent. 

We have also examined 8.4GHz A configuration obser- 
vations made by other observers. From an observation on 
1991 June 24, Ku kula et all (|l995h list the core peak flux as 
0.54 mJy and the integral flux as 0.60 mJy from a map with 
rms noise level of 73 /iJy/beam. In an observation on 1991 
September 01, which we analysed ourselves, we find a peak 
flux of 0.50 mJy with rms noise of 40 /xJy/beam (we do not 
quote the integral flux density as the beam is more distorted 
than that of any other A configuration maps which we con- 
sider here). These measurements are very similar to those 
which we have measured in our own later observations. 

We conclude that although we formally measure varia- 
tions of ~ 25 per cent on a timescale of 2 weeks between the 
first two observations in the maps made with fixed restor- 
ing beam, the variation is barely more than 0.1 mJy and, 
given the possible 30% larger error, the observations could 
be consistent with a constant source. 
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Figure 2. Top panel: A configuration peak 8.4GHz core flux 
densities with fixed beam (filled circles) and default/not fixed 
beam (open pentagons). Note the non fixed beam points (open 
hexagons) were moved forwards by 2 days to the make them eas- 
ier to see. The flux values were derived from jmfit by fitting 
to maps made with both fixed and default beams. See text for 
more details. 2nd panel: A configuration integral intensity val- 
ues for 8.4GHz core derived from maps made with a non fixed 
restoring beam. 3rd panel: RXTE fluxes slightly smoothed (run- 
ning gaussian with a half-width 2 days) to remove high frequency 
variability. Bottom panel: 2-10 keV RXTE fluxes (xlO -11 erg 
cm -2 s _1 ). Observations occur approximately once every 2 days 
with typical duration lksec. 



Table 2. A configuration core peak flux density values at 8.4 
GHz and 4.8 GHz for maps made with fixed restoring beams. 



Radio 


Peak Flux Density 


Peak Flux Density 


Date 


at 8.4 GHz 


at 4.8 GHz 


JD-2450000 


xlO- 4 Jy 


xlO" 4 Jy 


1852.222 


4.79 ±0.30 


5.91 ±0.36 


1865.185 


5.97 ±0.36 


6.36 ±0.37 


1879.148 


5.93 ±0.35 


6.59 ±0.38 


1893.131 


5.36 ±0.34 


6.41 ±0.39 


1909.066 


6.01 ±0.35 


7.39 ±0.43 


1929.926 


4.97 ±0.31 


5.98 ±0.36 



3.1.2 8.4GHz to 4.8GHz spectral variations 

We carried out the same map making procedure on the 
4.8GHz A configuration observations. In these cases we fixed 
the beam size at 0.385 x 0.328 arcsec with position angle -15 
degrees, i.e. quite similar to the PA of the 8.4GHz observa- 
tions. The resulting peak flux densities are listed in Table [2] 
and plotted in Fig [3] A broadly similar pattern of variations 
to that seen at 8.4GHz is found and the amplitudes of vari- 
ability are similar to those seen in NGC 5548 by [Wrobel 

In Fig[3]we also show the 8.4-4.8 GHz two-point spectral 
index, a, where S(v) oc v a . Given the difficulty of measuring 



core fluxes we do not over- interpret these results, however we 
do note that all the maps from which the fluxes were derived 
were constructed in an identical way. With that proviso we 
note that these observati ons reveal an averag e a = —0.2, 
similar to that found by Chris topoulou et al.l (1997), with 
a slight hardening at the beginning of the observations fol- 
lowed by a gradual softening such as might be explained by 
an injection of electrons with a relatively flat energy dis- 
tribution, e.g. perhaps as a result of a shock, followed by 
radiative energy losses. 

As this paper is concerned mainly with the radio vari- 
ability and as the problems of contamination of nuclear flux 
by extended emission are even greater at 4.8 GHz than at 
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Figure 3. Upper panel: The core peak flux density (mJy) at 4.8 
GHz and 8.4 GHz during the A configuration derived from maps 
made with identical restoring beams at each frequency. Middle 
panel: The 2-point 8.4 to 4.8 GHz spectral index, a where S(u) oc 
u a . Bottom panel: Smoothed X-ray flux for the time period of 
the observations in A configuration. 

8.4 GHz, we do not consider the 4.8GHz from configurations 
other than A configuration in this paper. However in the A 
configuration the nucleus is still reasonably resolved at 4.8 
GHz and so we do present those data here. A full discussion 
of the spectral morphology of NGC 4051 will be presented 
in a later paper. 

3.1.3 Radio /X-ray Correlation 

In Fig. [2] we also show (bottom panel) the 2-10 keV X- 
ray fluxes derived from our monitoring programme with 
the Rossi X-ray Timing Explorer (RXTE) using the Pro- 
portional Counter Array (PCA). During the period of the 
A configuration the RXTE observations typically occurred 
every 2 days, with a duration of lksec. Description of the 
X-r ay reduction tec hniques can be found in other papers 
fe.g lM c Hardv et alJliooi ). It is well known that NGC 4051 
varies very rapidly, and that rapid variability can be seen 
in Fig. [2] Although there are slow long term trends in the 
average flux of the source, obser vations with h igher time res- 
olution by XMM-Newton (e.g lM c Hardv et alJl2004h reveal 
that variations by a factor of 2 occur not uncommonly on 
timescales of an hour. It is therefore relevant to consider how 
close in time the X-ray and radio observations should be in 
order to be useful for any study of correlated variability. 

The answer depends on what we think the physical 
mechanisms responsible for the X-ray and radio emission 
are. For example, in blazars where the variable emission 
in both the X-ray and lower frequency bands is believed 
to come from shocks in a relativist ic jet oriented towards 
the observer, time separations of a few hours or less are 
required (e.g. lM c Hardv et alJl2007h . Although the VLBI ra- 



dio morphology of three compact co-linear components is 
most easily explained by the presence of an unseen jet, the 
radio/X-ray ratio is much lower than in typical blazars and 
a similar relativistic jet-like origin for both emissions is un- 
likely. The currently most favoured, and energetically sim- 
plest, paradigm for the X-ray emission is Comptonisation of 
optical- UV disc photons by very hot thermal (T > 10 9 K) or 
non- thermal electrons in a corona above the accretion disc 
(e.g. lShapiro et al"1ll976h . One possibility for the radio emis- 
sion is that it is synchrotron emission from that same corona, 
although the size of the radio emitting region at the frequen- 
cies considered here would be typically 100 times larger than 
that of the X-ray emitting region (see lLaor fc Beharll2008l . 
for an extensive discussion of this model). In that case we 
would not expect to detect radio variability on the rapid 
timescales seen in the X-ray band and it would be more 
appropriate to look for correlations between the radio emis- 
sion and a longer-term averaged X-ray emission which may 
represent longer term accretion rate fuelling of the overall 
emission regions. 

In order to avoid pre-judging the issue, we have 
therefore investigated the relationship between both the 
smoothed, and unsmoothed, X-ray fluxes. (The X-ray 
lightcurve is only gently smoothed with a running gaus- 
sian function of half- width 2 days.) The smoothed X-ray 
lightcurve is shown in the second panel from the bottom in 
Fig. El 

In case the time difference between the radio and X- 
ray observations is important, we list the radio observation 
dates and the nearest RXTE observation dates in Table [3] 
We list both the directly observed, i.e. unsmoothed X-ray 
fluxes and the fluxes interpolated to the time of the radio 
observations. We note that the greatest separation is 0.59d 
but that three of the separations are within 3 hours. We note 
that, in fact, there is not a great deal of difference between 
the nearest observed and interpolated X-ray fluxes. 

In Fig. [4] we plot the 8.4GHz radio flux from the maps 
made with the same restoring beam against both the in- 
terpolated unsmoothed and interpolated smoothed X-ray 
fluxes. There is a weak correlation between the radio and 
X-ray fluxes although we repeat our caution that, with a 
slightly larger error, the radio fluxes could be consistent 
with being constant. Smoothing the X-ray flux does not im- 
prove the appearance of the relationship. In Figs. [5] and [6] 
we plot the same relationship for the peak, and integral, 
8.4GHz flux densities derived from the maps made with the 
default restoring beams. The same general patterns are ob- 
served. We are able to measure simple power law X-ray spec- 
tral indices for each of the RXTE observations and, in Fig. [3 
we plot the interpolated photon number indices, from both 
smoothed and unsmoothed light curves as before, against ra- 
dio flux density for the 8.4GHz maps made with the same 
restoring beam. Again a similar pattern is seen although we 
caution that there is usually a good relationship between 
X-ray flux and spectral index for Seyfert galaxies such as 
NGC 4051 (jLamer et alJl20Q3h so these plots are not com- 
pletely independent of the previous flux-flux plots. 

Hardening of radio spectral with increasing flux followed 
by softening with decreasing flux can be simply explained in 
the context of jet models by the injection of particles with a 
hard spectrum followed by subsequent radiative losses, and 
so might explain our observations of radio variability. How- 
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Table 3. 8.4 GHz observation dates with the nearest RXTE X- 
ray observation dates. X-ray fluxes are given both as observed 
and interpolated to the time of the radio observation. Dates are 
JD-2450000. X-ray fluxes are in units of 10 -11 erg cm -2 s _1 . 



Radio 
Date 

1852.222 

1865.185 
1879.148 
1893.131 
1909.066 



Interpolated 
X-ray flux 

1.10 ±0.13 

1.84 ±0.14 
1.61 ±0.11 
1.23 ±0.10 
1.14 ±0.15 



X-ray 
Date 

1851.631 
1853.551 
1865.178 
1879.632 
1893.637 
1909.184 



Observed 
X-ray flux 

1.37 ±0.13 
0.51 ±0.11 
1.84 ±0.14 
1.83 ±0.13 
1.41 ±0.13 
1.14 ±0.19 






1929.926 0.35 ±0.17 1929.794 0.35 ±0.17 



^ 0.6 
^0.55 
3 0.5 
3 0.45 

. ~, i : : 

I n terp o II a i e d u n s n loothed X — r a y flux (erg c m - 



>, 0.6 
^0.55 

oo 

0.45 

0.5 1 : . 2 2.5 

Interpolated smoothed X-ray flux (erg cm -2 s -1 x 1 -11 ) 

Figure 4. 8.4GHz A configuration peak flux densities derived 
from maps made with the same restoring beam plotted against 
the X-ray flux. In the top panel we use the observed X-ray fluxes 
from the bottom panel of Fig. [2] interpolated to the time of the 
radio observations. The exact time of the radio and nearest X-ray 
observations are given in Table 2. All radio observations have an 
X-ray observation within 0.6d and, in 3 cases, within 3 hours. In 
the bottom panel here we plot the X-ray fluxes interpolated from 
the slightly smoothed X-ray fluxes (third panel of Fig. [2}. 
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Figure 5. Peak 8.4GHz core flux densities derived from maps 
made with the default restoring beam in A configuration plot- 
ted against the X-ray flux. In the top panel we use the observed 
unsmoothed X-ray fluxes, and in the bottom panel we use the 
smoothed X-ray fluxes. 
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Figure 6. Integral Intensity radio flux densities at 8.4GHz core 
flux derived from maps made with the default restoring beam in 
A configuration plotted against the observed X-ray fluxes (top 
panel) and the smoothed X-ray fluxes (bottom panel) . 



ever if the X-ray emission also arose predominantly in the 
jet we would expect a similar spectral behaviour, which is 
not what we see here. We refrain from speculating too much 
given the quality of the data but the correlated variability 
seen here, if real, may simply reflect changes in fuelling rate 
to both emission regions. 

3.2 During period of B configuration observations 

3.2.1 8.4GHz observations 

We see, from Fig.[T] that in the B configuration observations 
the core is not resolved from the two close components on 
either side. Measurement of the core flux therefore becomes 
more difficult and flux density measurement errors will be 
larger than for the A configuration observations. Nonetheless 
it is still possible to measure the peak, and integrated, flux 
densities. 

In Fig [8] we show the lightcurves during the period of 
the B configuration observations. Here the radio flux den- 
sities used are both the peak core flux densities measured 
from maps made with the default beamshape and those with 
a fixed beam . As with the A configuration observations, 
the values derived from maps made with the same/fixed 




0.45 0.5 0.55 0.6 

8.4GHz flux density (rnJy) 

2.5 I 1 1 1 1 1 1 1 n 1 ' 1 1 1 1 1 1 1 



0.45 0.5 0.55 0.6 

8.4GHz flux density (rnJy) 

Figure 7. Peak 8.4GHz flux densities from A configuration maps 
with the same restoring beam plotted against RXTE 2-10 keV 
photon number spectral index, T. In the top panel we interpolate 
from a slightly smoothed lightcurve of F and in the bottom panel 
we simply interpolate T between the nearest observed values. 
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Figure 8. In the top panel we plot the peak core flux densi- 
ties at 8.4GHz derived from B configuration maps made with the 
default /unfixed restoring beam (open pentagons) and fixed beam 
(filled circles). Note, the unfixed beam points have been increased 
by 2 days, to stop overlap and make them easier to see. We plot 
these against the unsmoothed and smoothed X-ray flux values or 
the period of the B configuration observations in the second and 
third panel respectively. 



restoring beam (i.e. 0.76 x 0.64 arcsec and position angle 
of -25.13 degrees) are quite similar to those made with the 
default beam. The X-ray observations are as described for 
the A configuration observations. We again see, at most, 
only marginal radio variability with the 8.4 GHz flux den- 
sity consistent with being constant. 
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Figure 9. Peak radio flux densities of the core at 8.4GHz derived 
from maps made with the default restoring beam in B configu- 
ration plotted against the observed X-ray fluxes (top panel) and 
the smoothed X-ray fluxes (bottom panel). 



array configurations there is no evidence of strong variability 
although, over months timescales, it appears visually as if 
there is a strong correlation between the radio and X-ray 
fluxes. 

Unfortunately, by chance, the period of the lowest X- 
ray flux occurs during the A configuration observations and 
the period of largest X-ray flux occurs during the D configu- 
ration observations. However note that here we are plotting 
the observed radio flux densities and have not yet removed 
contributions to the core flux density from extended emis- 
sion. As the nuclear radio flux will be most contaminated by 
extended emission during the D configuration observations, 
it is probable that at least part of the apparent correlation 
between the radio and X-ray fluxes is actually a consequence 
of the changing array configuration. In order to properly 
compare the radio and X-ray variability we must therefore 
remove the extended radio flux, as we describe below. 



3.2.2 Radio /X-ray correlation 

The relationship between the radio and X-ray fluxes in the 
B configuration is shown in Fig. [9] again displaying both 
the smoothed and unsmoothed X-ray fluxes. The hint of a 
positive X-ray /radio correlation seen in the A configuration 
observations is not seen here. With the possible exception of 
the highest X-ray flux point, there is no sign of any variation 
of the radio flux despite considerable variation in the X-ray 
flux. Unfortunately that highest X-ray flux point is derived 
from the largest extrapolation of fluxes as the RXTE obser- 
vation closest to that time had zero useful exposure. Thus 
we have had to extrapolate between two observations which 
were 4d apart. 

3.3 All configurations 

Although the errors increase, it is possible to measure peak 
and integrated flux densities also for the C and D array 
configurations. In FigHUJwe plot the resultant 8.4 GHz fluxes 
densities for all configurations together with the observed, 
and slightly smoothed, X-ray fluxes. Within the individual 



4 COMPARING RADIO FLUXES BETWEEN 
ARRAYS 

The procedure for removing the contribution to the core flux 
density from extended emission is carried out in stages, off- 
setting by one array configuration change at a time. For ex- 
ample, in offsetting from B configuration to A configuration 
we first determine, from real observations with B configu- 
ration, what the typical beam shape and size is in B con- 
figuration. We then take the A configuration observations 
and taper the uv datasets, using uvrange within IMAGR, 
so that the beam sizes are the same as they would be in 
B configuration. We refer to the resulting map as a pseudo 
B configuration map. We then measure both the peak and 
integral flux densities in the pseudo B array map and de- 
termine the differences, or offsets, between values measured 
from the original, untapered, A configuration map. We re- 
peat the process for all of the A configuration datasets and 
determine the average offsets to the B configuration, and an 
rms error thereon, for the peak and integral flux densities. 

To offset the real observed B configuration observations 
to an A configuration observation we remove the offset de- 
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Figure 10. Top panel - directly observed X-ray fluxes. Second 
panel - X-ray flux smoothed with a 4 point boxcar. Third panel 
-integral flux densities at 8.4GHz peak labelled with VLA ar- 
ray configuration. Fourth (Bottom) panel - peak flux densities 
at 8.4GHz peak labelled with VLA array configuration. The ap- 
parent variability in the radio flux densities is largely due to the 
changing array configuration of the VLA every ~ 100 days. 



Figure 11. Integral flux density (top panel) and peak flux den- 
sity (bottom panel) at 8.4GHz of the core. Flux densities for all 
observations are offset to the A array configuration. The peak flux 
density is consistent with being constant. 

flux densities. In the next section we simulate the offsetting 
method for integral flux densities to determine its accuracy. 



rived above from each of the observed B configuration flux 
densities and we combine the rms error on the offset with 
the observational error on the B configuration flux densities 
in quadrature. 

We perform this offset method on the other configura- 
tions to get average offset values for D to C, and C to B. 
Combining these offsets with the offset from B to A allows 
us to offset all data to an A array configuration. For ex- 
ample, the average peak flux density offset in going from A 
to B configuration is 0.32±0.07mJy and in going from C to 
D configuration it is 2.31±0.21mJy. In offsetting over more 
than one array configuration change we combine the offsets 
linearly and the errors in quadrature. Thus these average off- 
set values were used to offset all flux density values to those 
of an equivalent A configuration observation. The resultant 
offset peak and integral flux density lightcurves are shown 
in Fig. im In the intermediate configurations such as BnC 
there are few epochs of observation, and the configuration is 
slightly different in all of them. It is therefore very difficult 
to estimate the offset. These observations give rise to the 
more deviant peak flux densities. However even including 
these points the peak flux densities are consistent with be- 
ing constant as long as we assume the more realistic larger 
error. For the integral flux densities, although the variations 
are now considerably reducted, there is still some residual 
variability, which still follows the X-ray variability. However 
we note that the changes still mainly follow the changing ar- 
ray configurations. Within all configurations apart from the 
A configuration, there is no evidence for variability so one 
must suspect an error in the offsetting procedure for integral 



5 SIMULATING THE OFFSET METHOD 

The first step in simulating the offset method was to make 
a model of the source. Although it is preferable that the 
model resembles the real source, it is not essential that it is 
identical to the real source as the aim of the simulations is 
only to test the offsetting method and could, in principle, 
be carried out on any extended source. 

To make the model we combined all the uv data sets 
from the A, B, C & D configurations into one uv data set 
which we imaged (Fig. I12[) . The 28 clean components from 
this combined image which best represented the structure 
of the AGN were then selected. Using uvsub we added the 
model contribution from these components to the residual 
uv data sets produced by IMAGR during our original imaging 
of the 22 observations in each of the 4 main array configura- 
tions. The resultant data sets represent how our actual ob- 
servational setups would have seen our model of NGC 4051. 
We then repeated the offsetting procedure, which we de- 
scribed above for the real data, on the simulated data, to 
produce a simulated integrated flux density offset to the A 
configuration. This procedure should result in a constant 
flux density level equal to the model A configuration values. 
The actual simulated offset lightcurve is shown in Fig. 1131 
The dotted horizontal line in Fig. [13] represents the expected 
A array integral flux density value. We see that the B con- 
figuration data are at the same level as the A configuration 
data but the C and D configuration data are too high by 
about 0.5-0.6 mJy which therefore represents an approxi- 
mate systematic error in our process of offsetting the integral 
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Figure 12. Image made by combining the UV data from all 
single arrays combined at 8.4GHz and restored with a typical B 
configuration beam (0.65" x 0.6"). The rms noise level is 1.103 x 
10 -05 Jy/beam and contours are at the following multiples of the 
noise level -2, 2, 3, 4, 5, 6, 7, 8, 9, 12, 15, 20, 25, 30, 35, 40, 45, 
50, 55, 60. 

C and D configuration flux densities. If this systematic error 
is subtracted from the original offset integral flux densities 
(Fig. Hip , then most of the apparent long term variation is 
removed from the integral flux densities. It is unfortunate 
that the array configuration changes follow the approximate 
changes in long term X-ray flux, and hence that the very 
small residual variation in the integral flux densities appears 
to follow the X-ray flux. However given the particular diffi- 
culties of measuring the integral flux densities, the fact that 
the changes closely follow the array changes and the fact 
that long term variability is not seen in the peak flux densi- 
ties, we conclude that the apparent small residual long term 
variability in the integral flux densities does not represent 
real variation of the core flux. 



6 DISCUSSION 

6.1 Radio Variability 

We have monitored the radio and X-ray variability of the 
probable soft state Seyfert galaxy NGC 4051 over a period 
of 16 months in 2000 and 2001 with observations approxi- 
mately every 2 weeks. We have carefully subtracted the con- 
tributions to the core flux density from the differing amounts 
of extended flux density which will be detected in the dif- 
ferent restoring beams in the different VLA configurations. 
Thus we have produced lightcurves equal to those which 
would have been produced by observations all with the A 
array. Over that period the core peak flux density which in 
this context is the best measurement of true core radio lumi- 
nosity, is consistent with being constant. Also, the average 
peak radio flux density during our observations is, within 
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Figure 13. Integral flux densities for the model NGC 4051 data, 
offset to A array. Here the mean A configuration model value is 0.9 
mjy (shown by the dotted line). In offsetting from array to array, 
offsetting errors are added in quadrature to the jmfit error and 
the 5% error on the flux. We note that the B configuration model 
data, when offset to A are at the expected A array flux values, 
but the mean C and D configuration data differ, on average, from 
the expected offset value by ~ 0.5 mJy. 



the errors, the same as during two previous observations in 
June and September of 1991. 

Within the A configuration, which is the only config- 
uration within which the core is resolved from close neigh- 
bouring components to the east and west, there is a hint 
of correlated X-ray and radio variability, with radio (and 
X-ray) spectral variations approximately following the flux 
variations. However although the X-ray flux varies by fac- 
tors of a few, the maximum radio variations are, at most, 25 
percent between observations separated by ^ 2 weeks. How- 
ever the amplitude of variability is, at most, 0.12mJy. Whilst 
the percentage changes, timescales of change, and spectral 
hardening during a rise are entirely consistent with the 
observations of radi o variab i lity o f another similar Seyfert 
galaxy, NGC 5548 (|Wrobell liooch . NGC 5548 is approxi- 
mately ten times brighter than NGC 4051 and so there is 
no doubt about the credibility of its radio variability. Un- 
fortunately for NGC 4051 the credibility of the variability 
depends strongly on the error associated with each measure- 
ment. Taking the larger of the errors which could be asso- 
ciated with the radio measurements (discussed earlier) the 
radio lightcurve would be consistent with being constant. 
Thus overall we conclude that there is no strong evidence 
for any large amplitude variability in the core radio flux 
of NGC 4051 at 8.4GHz over the course of a year but we 
cannot absolutely rule out very low amplitude (~ 0.1m Jy) 
variability correlated with the X-ray variability. 
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Figure 14. The radio/X-ray /black hole mass 'fundamental 
plane' for hard state accreting black holes (i.e. h ard state X-ray 
binar y systems and liners) using the data from the lKording et ahl 
(2006) sample. NGC 4051 is the large black star, and is labelled. 
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Figure 15. Our mean values of X-ray and peak radio luminos- 
ity for NGC 4051 (l arge 5-pointed s tar) s uperposed on the lower 
panel from Fig. 2 of lLaor L Beharl (2008). The solid line is the 
Gndel -Benz relation, i.e. Lr/Lx — 10 -5 , for coronally active 
stars dGudel fc Bendll993h . 

6.2 Implications for emission models 

In Fig. [TJ] we plot our values of the radio luminosity and 
the average X-ray luminosity on the so-calle d ' fundamen- 
tal p l ane' of black hole accretion activity (Mer loni et al.l 
120031 : iFalcke et all 120041 ). Here we plot a version of the 
plan e using only 'hard state' black holes from the sample 
of Kording et al. (2006). We note that NGC 4051, which is 
almost certainly a 'soft state' object, lies about a decade 
below (in the sense of being more radio quiet) the best-fit 
relationship denned by the hard state objects, although it 
is just about within the scatter of the relationship. How- 
ever hard state objects typically show much more radio 
variability than NGC 4051 and are typically described by 
Lr oc Ly where B~ 0.7 for variations in in dividual objects 
([Corbel et al.ll2003l : iGallo et alJbOOl l2006h . For NGC 4051 
P — 0, or possibly ft — 0.1 during the A configuration ob- 
servations. 



Bell et al.l (|2010h find a significant correlation between 
X-ray and radio emission in the LINER NGC 7213 with 
factors of ~ 2 or more variability in both bands. Given 
its very low accretion rate, NGC 7213 is thought to be a 
hard state system so correla ted rad i o/X-r ay variability is ex- 
pected. However although Wr obeli (|200Q[ ) has demonstrated 
radio variability in the Seyfert galaxy NGC 5548, which has 
an accretion rate not far below that of NGC 4051 and so is 
probably also a soft state system, there have been no pre- 
vious studies of correlated X-ray /radio variability in simi- 
lar Seyfert galaxies. The X-ray/radio relations hip has, how- 
ever, b een investigated for samples of objects. IPanessa et al.l 
(|2007T ) find p = 0.97 at 4.8GHz and f3 = 0.98 - 1.25 at 
15GHz for a sample of Seyfert galaxies. For a sample of 
ROS AT-selected AGN and using 1.4GHz radio observa- 
tions, |lHnkmanneian([2QQ3) find f3 = 0.48±0.049 for radio 
quiet objects, similar to Seyferts, (and f3 = 1.012 ±0.083 for 
radio loud objects). Thus although there is a wide range in 
the reported values of /?, there is good evidence for radio/X- 
ray correlation within samples of Seyferts and radio quiet 
AGN. However within the only well studied 'soft state' 
Seyfert so far, i.e. NGC 4051, there is no evidence for a 
value of f3 much above 0.1. 

One possibility is that even in the VLA A configuration 
8.4GHz maps, the core still contains substantial contribu- 
tion from emission which is extended but on scales below 
that detectab le with the VLA. We not e that with a beam 
of ~ 8mas, iGiroletti Panessal (|2009l ) find an unresolved 
source with 5GHz flux density of 0.20 =b 0.02mJy, which is 
substantially less than our average 4.8GHz measurement of 
no 0.6mJy. Our tentative observations of ~ 0.1m Jy variabil- 
ity would then provide a value of f3 consistent with the val- 
ues derived from observations of samples of similar sources, 
although NGC 4051 would then be even more radio quiet 
compared to the hard state objects which define the funda- 
mental plane. 

The VLBI observations show structure which is strongly 
indicative of a collimated jet ([Giroletti &; Panessal [2009; 
lM c Hardvll2005h . (In a future paper we will present our 
higher resolution global VLBI observations which confirm 
this structure.) Thus as first sight the lack of radio vari- 
ability is surprising. However if the jet is in the plane of 
the sky, rather than pointing at t he observer as in b lazars, 
radio variability would be slower. Falcke et al . (2001) spec- 
ulate similarly regarding the lack of strong radio variability 
in radio-loud quasars although they do not reach any firm 
conclusions. 

Another alternative is that the radio emission in Seyfert 
galaxies does not come from a jet, as in hard state 
sources, but arises in the outer parts of the X-ray emit- 
ting corona (|Laor &; Beharl 120081). In Fig. [T5l we plot our 
current data onto Fig. 2 from lLaor &; Beharl |2008) an d see 
that NGC 4051 agrees reasonably with the Giidel-Benz ll993l 
relationship for coronally active stars, i.e. Lr/Lx = 10 -5 . 
In the coronal model, rapid large amplitude radio variability 
is not expected at frequencies where the source is optically 
thick. This situation is approximately the case during our A 
configuration observations which indicate a spectral index 
of a no —0.2. Such variability is, indeed, not seen. Similar 
considerations would, of course, also apply to some extent 
to a side-on jet and the highly collimated co-linear structure 
of the compact VLBI components is a significant problem 
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for the coronal model. As jets are generally agreed to arise 
from the corona, the truth may lie part way between the 
pure coronal and pure jet models. In this galaxy where an 
outflow of large solid ang le is seen in [OIII] observations 
(|Christopoulou et al.l ll997). the compact radio components 
separated from the nucleus may arise from a central, higher 
velocity region within the overall outflow. The core radio 
emission may arise from a combination of a jet and a corona. 
If there is any long term correlated variations between the 
radio and X-ray bands, they may simply reflect long term 
changes in accretion rate, upon which both bands depend, 
rather than a direct link between the emission processes. 



7 CONCLUSIONS 

Despite factors of 10 X-ray variability, we see no clear 
evidence for variability of the core of the soft state narrow 
line Seyfert 1 galaxy NGC 4051 at 8.4GHz over a 16 
month period of monitoring at 2- weekly intervals, with 
the possible exception of very low amplitude (~ 0.12mJy) 
variations during the A configuration observations where 
the core is best resolved from surrounding structures. The 
latter tentative variations correlate weakly with the much 
larger amplitude X-ray variations. Our resultant radio and 
mean X-ray luminosity make NGC 4051 about a decade 
radio quieter than the hard state objects which define the 
'fundamental plane' for hard state accreting black holes, 
although the scatter about the plane is almost of the same 
order. Given the collimated VLBI structure which hints at 
the presence of an unseen jet, the lack of radio variability 
is, at first sight, surprising, although a side-on jet, observed 
at a frequency where it is optically thick, would not vary 
as rapidly as a face-on blazar jet. A coronal model agrees 
well with the radio/X-ray flux ratio, and the lack of radio 
variability, but the collimated radio structure is then hard 
to explain. A mixture of jet and coronal emission may 
explain the observations but further, more sensitive, radio 
observations, with high angular resolution and a fixed 
beamshape, are required to confirm whether the radio 
emission from NGC 4051 does indeed vary. 
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